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Abstract In this work we consider how the ap¬ 
pearance of gradients of pseudoscalar condensates 
in dense systems may possibly influence the trans¬ 
port properties of photons in such a medium as 
well as other thermodynamic characteristics. We 
adopt the hypothesis that in regions where the 
pseudoscalar density gradient is large the proper¬ 
ties of photons and fermions are governed by the 
usual lagrangian extended with a Chern-Simons 
interaction for photons and a constant axial field 
for fermions. We find that these new pieces in the 
lagrangian produce non-trivial reflection coeffi¬ 
cients both for photons and fermions when enter¬ 
ing or leaving a region where the pseudoscalar has 
a non-zero gradient. A varying pseudoscalar den¬ 
sity may also lead to instability of some fermion 
and boson modes and modify some properties of 
the Fermi sea. We speculate that some of these 
modifications could influence the cooling rate of 
stellar matter (for instance in compact stars) and 
have other observable consequences. While quan¬ 
titative results may depend on precise astrophysi- 
cal details most of the consequences are quite uni¬ 
versal and consideration should be given to this 
possibility. 
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1 Introduction 

Interest in possible violations of Lorentz and other 
fundamental symmetries emerged P) after the pa¬ 
per [2], where electrodynamics with an additional 
Chern-Simons (CS) parity-odd term including a 
constant CS four-vector was considered. The same 
authors entertained the possibility that such an 
effect could already be visible in the radiation ob¬ 
served from distant radio galaxies, arriving at a 
negative conclusion. Nevertheless, there are some 
areas where modified Maxwell-Chern-Simons (MCS) 
electrodynamics (also known as Carroll-Field-Jackiw 
electrodynamics) may be relevant. 

It has been suggested that parity breaking 
phenomen^ may take place in peripheral heavy 
ion collisions [3] manifesting itself in the so-called 
Chiral Magnetic Effect (CME)[4]. Spontaneous 
parity violation might also occur for sufficiently 
large values of the baryon density [5]. Recently 
several experiments in heavy ion collisions have 
also indicated an abnormal yield of lepton pairs [6l 
[71 and it was conjectured that the effect may be 
another manifestation of local parity breaking in 
colliding nuclei [H[9|. This effect would be due to 
an interplay between topological fluctuations in 
QCD (see [TOl[TT] and refs, therein) and MCS elec¬ 
trodynamics. 

Another source of possible macroscopic mani¬ 
festations of parity breaking is the presence of an 
axion background. Even though there is no di¬ 
rect evidence of the existence of axions yet it has 
been speculated that parity odd regions may oc- 

1 Parity breaking in strong interactions is possible in 
out-of-equilibrium processes or with a non-zero chem¬ 
ical potential. 
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cur after Bose-Einstein condensation of axion or 
axion-like fields m-m at several astrophysical 
scales. In particular this could take place in the 
interior of stars |16j. Recent investigations specu¬ 
late that these axion condensate^ could be rather 
compact in size m- 

In dense nuclear matter one should contem¬ 
plate also phenomena such as neutral pion conden¬ 
sation]]]))] or the occurrence of a disoriented chiral 
condensate [20:. Both phenomena have received 
considerable attention in the past. 

We will generically refer to these phenomena 
as pseudoscalar condensation and most of its con¬ 
sequences will be independent of the precise mech¬ 
anism triggering the appearance of local parity 
breaking. 

In fact what matters for our purposes is not 
the existence of pseudoscalar condensates per se 
but rather gradients of the pseudoscalar density. 
A region with a strictly constant pseudoscalar 
density has a vanishing CS four-vector and its 
electromagnetic properties are described by usual 
Maxwell electrodynamics. On the contrary in re¬ 
gions where the pseudoscalar density is space- 
dependent a non-zero CS vector appears. In this 
region the relevant lagrangian to describe elec¬ 
tromagnetic interactions is the MCS lagrangian. 
Similar arguments can be applied to the constant 
axial-vector field coupling to fermions [21]. We shall 
therefore be interested in this case. 

In the next sections we discuss several effects 
resulting from the appearance of pseudoscalar gra¬ 
dients in dense systems. We first consider how 
they may affect the propagating modes of pho¬ 
tons and fermions. For the former, different polar¬ 
izations are affected differently, while in the case 
of fermions the Fermi sea is split in two, with 
slightly different Fermi levels. We then proceed 
to determine transmission and reflection coeffi¬ 
cients for photons crossing a layer where a pseu¬ 
doscalar condensate has a gradient, concluding 
that it may have a relevant influence in process 
mediated by radiative cooling if such gradients 
are sufficiently large compared to the tempera¬ 
ture. A similar phenomenon occurs in principle 
for fermions, but it is much suppressed due to 
their larger mass. Fermions can ‘decay’ (i.e. move 
from the higher Fermi sea to the lower) in such 
a medium emitting a photon mm- Conversely, 
photons propagating in such a medium but free 

2 These axion condensates should not be confused with 
the classical cold axion background — a viable candi¬ 
date for dark matter [lSj . 


otherwise can materialize in a fermion-antifermion 
pair in certain circumstances. To conclude we spec¬ 
ulate about the possible relevance of the above 
phenomena in degenerate fermion systems such as 
compact stars where, as previously argued pseu¬ 
doscalars condensates may be present. 

2 Spectrum of photons in a pseudoscalar 
condensate gradient 

Let us commence our analysis by assuming that 
in a dense system there are domains with differ¬ 
ent pseudoscalar densities. These domains may 
be either of axion or pion-like type. We will char¬ 
acterize the different domains by approximately 
constant values of the condensate a*. Obviously 
a = 0, i.e. no parity breaking at all, is also a pos¬ 
sibility. Typically, transition regions between dif¬ 
ferent domains and between a parity breaking re¬ 
gion and the normal vacuum will be present. The 
figure shows the generic situation we are consid¬ 
ering: region 2 is a transition region with non-zero 
pseudoscalar gradient separating regions 1 and 3. 
Inside a transition region —one with a non-zero 



Fig. 1 A sketch of possible domains with pseudoscalar 
condensates. 

gradient— the behavior of the photons can be de¬ 
scribed by Carroll-Field-Jackiw model, or MCS 
electrodynamics, with the Lagrange density 

C = - \ F a P(x)F a p(x) - i F^(x)F^(x) a(x) 

+ \ A^(x)m 2 t V^A^x), (1) 

where A^{x) and a(x) stand for the vector and 
effective background pseudoscalar fields respec¬ 
tively, F^ = \ e F P c is a dual field strength 
and rat is an effective transverse photon mass [231 
[24] (see Sec.6) which is generated by gauge invari¬ 
ant polarization tensor in plasma, containing the 
transversal projector = rf‘ v — d^d v /d 2 . If one 

uses the low momentum limit k —> 0, at high tem¬ 
peratures T and chemical potentials /i, this mass 
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happens to be constant, m 2 ~ ^T 2 + ^ 

0. 

The classical background is different for dif¬ 
ferent regions. In the area outside the boundary 
layers (regions 1 and 3) a(x) is constant. In this 
case the second term in Lagrangian does not give 
any corrections to field equation and one recov¬ 
ers standard electrodynamics. The only nontriv¬ 
ial area to describe is the region 2. We assume 
that in this layer a changes linearly from a_ to 
a + across the gap. Then inside the region 2 the 
relevant pseudoscalar background can be locally 
described by 

a(x) = C • x[0(C ■ (x - x _)) - 0« • {x - £+))] -> 
C-z = C^; (2) 

in which a fixed constant four vector with mass 
dimension is used as an argument. is actually 
proportional to the local gradient of the pseu¬ 
doscalar condensate. We assume that the wave 
function of photons is considerably shorter that 
the thickness of the layer and accordingly take the 
width of the latter to infinity to simplify our cal¬ 
culations. In this case Lorentz invariance would 
be violated in the Minkowski half space. 

If we consider a small area that does not feel 
the curvature of the shelH, we may choose for 
simplicity the first coordinate along the local ra¬ 
dius of curvature of the bubble and assume the 
first boundary X- to be located at x\ = 0. In this 
particular case (£ M ) = (0, 0, 0), 

d(x) = Czi 0(xi), (3) 

where we have now assumed that the thickness of 
the boundary is much larger than the character¬ 
istic photon wave length and accordingly we have 
taken x+ —> oo. This approximation makes the 
calculation simpler as it allows to decouple the 
effect of the two successive interfaces. 

Then the wave equation reads 

(□ + m 2 t )A v + Ce luap 0{x 1 )d a A p = 0 (4) 

Its solution in the MCS medium can be found by 
proper projection on the longitudinal and chiral 
modes (see | Appendix A[ ). The corresponding dis¬ 
persion relations are 

kiL = fc° = Vv 2 - k'j_; _ 

= \] w 2 - ml - k\ - CiV ^ 2 ^ ( 5 ) 

= \J uj 2 — m 2 — + (iy/oj 2 — k‘j_ , 

3 Again, this is a good approximation for photons 
whose wavelength is much shorter than the character¬ 
istic sizes involved 


where the index ‘0’ labels the medium with the 
usual vacuum dispersion relation, w is a photon 
energy and kj^ = (0, Aq, Aq) is a transversal pho¬ 
ton wave vector. This expression describes three 
different physical polarizations in the MCS medium 
for x\ > 0. Of course the standard Maxwellian 
behavior is recovered by setting £ = 0 in all ex¬ 
pressions, which corresponds to the region aq < 0, 
where there is no Lorentz symmetry breaking. 


3 Spectrum of fermions in a pseudoscalar 
condensate gradient 

In this section we outline the key points of the 
properties in 3+1 Minkowski space-time of free 
spinor fields in the presence of a Lorentz covari¬ 
ance breaking kinetic term associated with a con¬ 
stant axial-vector This vector is supposed to 

be generated in the region between two bubbles 
with different pseudoscalar condensate and its pres¬ 
ence will change drastically some fermion proper¬ 
ties in the vicinity of Fermi surface. 

The free fermion spectrum can be obtained 
from a modified Dirac equation in the momentum 
representation 

- m - 7^75) 4 ’ = 0 . ( 6 ) 

The solution of this equation is discussed 
in |Appendix B| The fermion spectrum is deter¬ 
mined by the on-shell condition 

(p 2 + b 2 - to 2 ) 2 + 46 2 to 2 - 4(6 • p) 2 = 0. (7) 

This equation has real solutions for any value 
of b However, a consistent quantization of the 
spinor field can be performed if there are two pairs 
of opposite-sign roots of eq. © and a mass gap 
between them. This condition holds true [25] for 
sufficiently small and a mass gap between posi¬ 
tive and negative frequencies exists provided that 
no solutions with po = 0 exist. Such solutions 
never arise for space-like b ^ which we employ in 
this paper. Moreover the choice b 2 < m 2 is quite 
plausible as the magnitude of the constant axial 
vector associated to a gradient of pseudoscalar 
condensate is expected to be of the order of a few 
keV at most (see discussion in next section). 

For the purely space-like case one can fix b M = 
(0, b) = (0, b, 0,0 ):b > 0 by the proper choice of 
coordinate systenQ Then the dispersion law is 

4 In concordant frames the results must be modified 
but the physical consequences remain the same 1261 
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given by 

uj 2 = p±+ Pi+b 2 + m 2 ± 2b\Jp\ + m 2 

= pi + (b ± \jpl + m 2 ^j ; 

P = (Pl,P2,P3); P± = (0,P2,P3)- (8) 


These solutions are separated by the stability cone. 
The stability border p 2 = 0 is described by 


\pi 



(9) 


4 Fermi sea in a pseudoscalar condensate 
gradient 


In this section we shall assume that the momenta 
involved fulfill the condition p -C in and therefore 
a non-relativistic approximation is valid. 

Taking this into account the dispersion re¬ 
lations for fermions in a medium with a non¬ 
vanishing pseudoscalar gradient would approxi¬ 
mately read 

u = m±b + & + EL. ( 10 ) 

2m 2m 

Let us assume that the medium is in thermal equi¬ 
librium with a temperature T. The two different 
mass-shell conditions established in the previous 
section will give rise to two different Fermi seas. 
These in the non-relativistic limit correspond to 
the ± branches in the above equation. Assuming 
that initially N+=N-=N the corresponding chem¬ 
ical potentials p± are obtained by solving eqs. , 


N =V duj 


\J (2m) 3 (uj ± b) 
4-7T 2 


exp —+ 1 


(ii) 


in respect to p. Clearly this gives p+ = p- + 2b 
and therefore ej > ej, (we assume that b is posi¬ 
tive; otherwise the situation is reversed). The ex¬ 
istence of pseudoscalar regions with a non-zero 
gradient thus naturally lead to a splitting of the 
Fermi sea in two, with a difference in levels di¬ 
rectly proportional to the gradient itself. For sim¬ 
plicity we are implicitly assuming that the gra¬ 
dient is constant (i.e. the pseudoscalar density 
varies linearly along the T’ direction). 

At least theoretically one could contemplate a 
single pseudoscalar condensate occupying a large 
part of a star primarily made of degenerate fermions 
(such a neutron star or a white dwarf) in its in¬ 
terior. Let us also assume for the sake of the ar¬ 
gument that the condensate varies linearly with 
the radius. From the previous considerations this 


would mean that the degenerate electrons or neu¬ 
trons that sustain the star from collapsing, in 
fight against the gravitational potential, are split. 
Half of the fermions (the “+’ ones) have to deal 
with more than their fair share of the weight of 
the star. The natural tendency under the action 
of gravity would then be to equilibrate both Fermi 
seas. In order to do so, some of the “+’ fermions 
have to turn into ones. Is this possible? The an¬ 
swer is yes and it will be discussed in the coming 
sections. The conversion would necessarily imply 
a large emission of photons with energies up to 
uj ~ 2b. 

The same phenomenon can present itself in 
smaller regions of the star leading to local im¬ 
balances in the Fermi sea that could trigger small 
starquakes with their corresponding photon emis¬ 
sions. 


5 Propagation of photons across 
pseudoscalar gradients 

Now we describe the propagation of photons in¬ 
side every of three regions indicated in Fig. 1. 
Regions 1 and 3 correspond to taking ( = 0 while 
C 7^ 0 in the intermediate transition region 2 
where the pseudoscalar condensate has a gradi¬ 
ent. In this section we neglect an transverse pho¬ 
ton mass TOt. As it is manifest in the previous 
equations, different regions lead to different dis¬ 
persion relations. As a consequence non-trivial re¬ 
flection and transmission coefficients between the 
different regions appear. This issue was first dis¬ 
cussed in m- 


5.1 Entering the boundary layer 

We assume that the thickness of the intermedi¬ 
ate shell (region 2 in Fig. 1) is much larger than 
a typical wavelength and a mean free path. Un¬ 
der these assumptions I® works. Photons fall on 
the boundary from a region where ordinary elec¬ 
trodynamics holds and attempt to penetrate in 
one governed by MCS electrodynamics. Match¬ 
ing conditions for this problem were discussed in 
detail in [25\ To understand which photons pen¬ 
etrates into the shell it is worth to recall the dis¬ 
cussion. Solution of can be found by using 
Fourier transformation over all components but 
x± and may be written in the form 

A v = 


( 12 ) 
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u v ^.e iklX1 +u v ^e x\ < 0 ; 


-zfeV 


E 


A L 


,CS cs 

0 lKl a X 1 _L /Ti . A X 1 


v vA ^e 1A Xl +v vA ^e~ 


, x L > 0 , 


where v and u depend on (to, The first 

index of v denotes the corresponding component 
of A v , v = 0,2,3, the second index A stands for 
the polarizations L, +, — and the arrows — 
point out the direction of particle propagation. It 
is necessary to distinguish among the various po¬ 
larizations because they obey different dispersion 
relations. Note that we contemplate the possibil¬ 
ity of photons developing a transverse mass nit 
which can be easily reinstated in the formulae. 

Then we take the initial amplitude, and us¬ 
ing eqs. (IC.23IIC.23l) in |Appcnclix C we find the 
amplitudes of the transmitted waves with polar¬ 
izations L, +, — 


‘I ± = 


1 


1 + 


1 T W 


(13) 


Longitudinally polarized photons are not affected 
by the change in the medium. 



Fig. 2 The geometry of photon propagation, n is a 
normal vector. In the left region (region 1) the pseu¬ 
doscalar condensate takes a constant value. Region 2 
(right) is assumed to be governed by MCS electrody¬ 
namics describing a varying pseudoscalar condensate 
as befitting a transition region 


The direction of outgoing photons corresponds 
to the angle /3 (Fig. 2). After the propagation 
through the boundary, their direction will be 
changed (angle a ) accordingly to their polariza¬ 
tion. We decompose k = k n n + k±r (Fig. 2). 
From [28| we know that k± remains the same af¬ 
ter crossing the boundary. This means that the 
trajectory of the particle and normal vector lie in 
one plane. However, k n changes. We introduce the 
new variable k^ s to describe the normal compo¬ 
nent of k after crossing the boundary. For x\ < 0 
one has 

k° 

cos/3 = — (14) 

to 

We are interested in finding the overall flux of out¬ 
going photons from region 1 to region 2 in Fig. 2. 
To do this we first consider a small volume near 


the boundary. We assume that it radiates uni¬ 
formly in all directions. For us is important the 
flux of energy which propagates outwards so we 
are interested only in upper half-sphere in Fig. 3. 
We assume that this small volume radiates with 



Fig. 3 Radiation of a layer 


a certain frequency Nq per solid angle. At this 
moment we neglect the photon mass to get some 
quantitative results. In order to find the total lu¬ 
minosity we should integrate over the solid angle, 
frequency and surface of the layer. The last inte¬ 
gration would be the same for the case with or 
without pseudoscalar condensation 

Lex / cIloNq f 2 d^ + (/3,to,0 

JO J min(-^ ,arccos(^-)) 

+ r dtoN% r d/3T-(/3,cc,C) (15) 

Jo Jo 

Here one can see that integration over angles in 
the first term begins from the value cos/3 = £. 
This value comes from kinematical condition of 
the positive polarization. It is easy to see from 
(1131) . that for T + in denominator we have a neg¬ 
ative value under the root if £ > k\. Physically 
it means, that for falling photons with k\ < C, 
it is kinematically forbidden to convert into pos¬ 
itive polarization photons in the medium with a 
linearly varying pseudoscalar field. For negative 
polarization there is no restriction, in (j5j fci_ is 
positive for any values of £ and as a result, we do 
not see any special limits of integral for negative 
polarization. 

Next we assume that the medium is in ther¬ 
mal equilibrium with a temperature T and ac¬ 
cordingly Nq oc u 3 /(eT — 1). We will compare 
this value with the luminosity of the same volume 
without any parity breaking boundary effect. Let 
us call this last value Lq 

( 2 dp. (16) 

0 

Finally, we can plot a graph which demonstrate 
the effect of the intermediate shell on the lumi¬ 
nosity. From Fig. 4 one can see that the effect of 


L 0 oc 


duo Nr 
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T 

Fig. 4 The outgoing energy flow ratio for the two cases 
(with and without pseudoscalar condensate) 


changing the dispersion relation across the bound¬ 
ary between region 1 with a = constant and re¬ 
gion 2 where a depends linearly in X\ is very 
noticeable. At large values of £ (compared with 
temperature) most photons are reflected from the 
boundary and accordingly the energy radiated 
across the boundary decreases. Note that as a 
consequence the radiation emitted does not cor¬ 
respond exactly to that of a black body with tem¬ 
perature T. 


We are interested only in transversal polariza¬ 
tions (we deal with photons, however, the longi¬ 
tudinal one does not feel the boundary anyway), 
so in our terms we write, 


T ± = 


2 kSi 


2cot(a) 


( 22 ) 


+ k n cot a + cot(fl) 

For our purposes it is necessary to express T 2 * 2 as 
a function of /? that means expressing a in terms 
of (3. For constant a(x), ui = |k| and one can use 
k n = oo cos a. So, m gives 
a; cos a 


cot a = cot (3 


\J a; 2 cos 2 a =F cos a 


(23) 


Where as usual the stands for different polar¬ 
izations. 

Solving this equation one can find the expres¬ 
sion for cot /3 for different polarization and the 
value of “T 2 * 2 , 


3^03,C,w) = 


(24) 


5.2 Escaping from the boundary layer 


After escaping the first region photons appear, if 
they are not reflected, in the intermediate shell 
where MCS electrodynamics is at work. To leave 
this medium and gain access to another domain 
where a = constant (possibly zero) photons have 
to pass through one more boundary. This corre¬ 
sponds to the boundary between regions 2 and 3 
in Fig. 1. We use the same technique as in the 
previous section. Figs. 2 and 3 still apply but re¬ 
versing the areas where ordinary and MCS elec¬ 
trodynamics apply. In this case 
uCS u 

-f- = cot(a); = cot^). (17) 

Furthermore we know that for spatial CS vec¬ 
tor there are two transversal polarizations in the 
MCS medium with the dispersion relations 


k nl = \/w 2 - kl =F Ct /^ 2 - 


V _L 


(18) 


or, since we consider photons and k n = y/u 2 — k\, 
k%l = y/k* T Ckn, (19) 


and 

kn ,, S \/ k n T C kn , / m 

— = cot(a); - — -= cot (/3). (20) 

Using the results obtained in [2S], one may find 
the transmission coefficient of outgoing particles 
for every polarization 




2 *Sf 

JcCS I b.0 ' 
^1 A ^1 


( 21 ) 


2cot/3 

cot g + ±C + VC 2 +4 u ,2 CQt^ P( l + COt 2 ?) 

yj 4a; 2 (1+cot 2 /3)-2C 2 l=2f y/C 2 + 4 ui 2 cot 2 /3(l+cot 2 /3) 

Using this formula we may find for any angle 
/3 the fraction of incoming photons succeeding in 
crossing the boundary at x+, i.e. the boundary 
between the regions 2 and 3 in Fig. 1. Like in 
the previous section we consider a small volume 
which radiates at certain frequency Af£) in a unit 
of solid angle and write total luminosity 

L oc J^° cLuiNq 2 d/3T + 03,w,C) 

+ f° duN% [ 2 dFZ-(P,u,Q (25) 

J o Jo 

It is worth commenting on the integration region 
of both terms. Looking at the expression (ED one 
may see that for negative polarization there are 
no any restriction on the kinematics of photon. 
This fact means we should integrate over all en¬ 
ergies of outgoing photons. In case of positive po¬ 
larization ui cannot be less than £. 

Like in previous section, to show the qualita¬ 
tive effect of pseudoscalar condensate we assume 
that the medium has a temperature T. From Fig. 
5 one can see that for £ > T there is a clear 
effect on the flux of outgoing energy. However, 
the boundary between regions 2 and 3 (Fig. 1) 
does not play a crucial role in the decreasing of 
energy flux effect (compare Fig. 4 and Fig. 5). 
Nevertheless, to get the total effect, one should 
take into account the impact of both boundaries. 
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Fig. 5 The outgoing energy flow ratio for the two cases 
(with and without pseudoscalar condensate) 
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Fig. 6 The outgoing energy flow ratio calculated for 
both boundaries 


If we neglect the thermal capacity of the shell 
H then we may plot a graph (Fig. [6]) showing 
the combined effect of two boundaries of the do¬ 
main. One sees that the effect may be substan¬ 
tial and in this case the thermal evolution of the 
bulk would be affected: radiative cooling (assum¬ 
ing that region 1 is hotter than region 3) would 
proceed more slowly. In systems like neutron stars 
or white dwarfs radiative cooling is not the domi¬ 
nant cooling mechanism and therefore the impact 
of the previous effect is limited. 


6 Fermion decay in a pseudoscalar 
gradient 

Fermions crossing domains with a pseudoscalar 
gradient will reveal the same effects of partial re¬ 
flection as the photons show for exactly the same 
reasons. The discussion of these effects is given in 
|Appcndix D| However if 6 <C m then the effects 
are clearly subleading. 

A more important phenomenon is the instabil¬ 
ity of fermions in a parity breaking environment. 
Namely, a fermion (be it an electron or a neutron) 
with a larger effective transversal mass m+ = 
VpI + m 2 + b may well decay into a fermion with 
a smaller effective transversal mass 
TO_ = 'M + m 2 — b by emitting a photon (in 
the case of electrons) or a neutral virtual pion 
(in the case of neutrons) that would immediately 
decay into two photons. Since the formalism is 

5 This implies assuming that the thickness of the shell 
is much smaller linear sizes of a typical domain where 
a = constant. 


identical in both cases let us restrict ourselves 
to the case of photon emission via this mecha¬ 
nism for distorted electrons. In this case the emer¬ 
gent photon will have a transversal mass m± = 
\/kf + m 2 + C 2 /4 ± C/2 where m t is the effective 
transversal photon mass in the plasma (if any). 

Let us restrict ourselves to the purely space¬ 
like 6 M = (0,6,0,0) and take b > 0. Besides one 
can use the remaining small Lorentz invariance 
0(1,2) in the hyperplane (0,2,3) and select out 
the transversal rest frame for the decaying elec¬ 
tron. Then the decay is possible for an electron 
with transversal mass m+ and momentum pp into 
an electron with transversal mass m_ and mo¬ 
mentum qp and a photon with momentum k M and 
transversal mass m± . Energy-momentum bal¬ 
ance for the decay, p M , consists of the 

following relations 

b + y Jp\ + TO 2 = + {yj 9i + m 2 - b) 2 + 

+ \jk\ + (yj k\ + m 2 + e 2 C 2 /4 + eC/2) 2 ; 

Pi = qi + h; q_L = -kj_, (26) 


where e = 0, ±1 for photon longitudinal and two 
C polarizations. 

In order to estimate the threshold of the re¬ 
action that allows *+’ electrons to decay into 
electrons, we neglect the positive contributions of 
the transverse momenta in the right-hand side, 


b + y]p\ + to 2 > yj'ql + m 2 - b + 

+ yjk\ + m 2 + e 2 C 2 /4 + eC/ 2 , (27) 

it is assumed that m > b. Thus, 

2b + yjpi + m 2 - eC ,/2 > m + \J+ e 2 C 2 /4, (28) 

wherefrom one derives the necessary inequalities 
to allow the decay to a particularly polarized pho¬ 
ton 


£ = 0 

£ = +1 


£ = —1 


26 + \/p\ + m 2 — m > m t ; 


(26 + y/p\ + to 2 — m) x 
x. /1 — 


C 


2b+ yjp\+m 2 —m 

(26 + \/p\ + TO 2 — m) X 
1 + 


> m t ; 


2b+^/pf+m 2 — m 


> m t ; 


(29) 


Evidently the absolute bound for the decay thresh¬ 
old is reached in the rest frame of a fermion p\ = 

0. At very high temperatures T and densities (chem¬ 
ical potentials p) and moderate photon momenta 
the effective photon mass[251 m 2 + 
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suppresses the transformation of l +’ electrons into 
electrons in the parity breaking layer. 

The transformation from one type of fermions 
to the other via photon emission is the process 
that may allow chemical equilibrium between the 
respective Fermi levels, eventually leadin g to a 
small excess of fermions in the systen 0 . As 
we have seen the viability of the process and the 
characteristics of the photon emitted depends cru¬ 
cially on the relation between the three parame¬ 
ters b, £ and uit- If the transition turns out to be 
energetically impossible because mt is too large, 
in the presence of this phenomenon a star has to 
be supported by a combination of the two Fermi 
seas. Notice that while in vacuum all “+’ electrons 
would like to decay to electrons, this is only 
possible for a limited number of electrons in the 
Fermi sea, as the new would have nowhere to 
go all levels being occupied. Therefore the contri¬ 
bution of this anomalous process to star cooling is 
very small in absolute terms although the sudden 
bursts of photons that may produce could rep¬ 
resent a clear observational signal. Whether this 
is the case or not depends on the transparency 
of the crust to these photons. If transmission is 
negligible they would only contribute to a slight 
increase of the star temperature. It is to be ex¬ 
pected though that a sudden reordering of the 
Fermi levels could produce some sort of glitch if 
the neutron star turns out to be a pulsar. 


7 Photon decay 

Finally let us discuss another phenomena, which 
may give contribution to the flux of outgoing par¬ 
ticles, namely, the possibility of high-energy pho¬ 
ton decay in a volume where a ^ constant. For 
a star thermodynamic balance the importance of 
this phenomenon is probably irrelevant but it is 
interesting on its own nevertheless. 

It was shown in [29] that a photon of positive 
polarization may decay in presence of a gradient 
of the pseudoscalar field into an e + e— pair. This 
process will suppress the number of outgoing pho¬ 
tons with positive chirality and possibly increase 
the number of outgoing electrons, positrons and 
antineutrino due to the process 

e + + n —> p + + v e . (30) 

6 The possible backreaction of this effect on the value 
of the pseudoscalar condensate itself has not been con¬ 
sidered. 


We will present here some calculations to under¬ 
stand what may be the quantitative effect of the 
photon decay. The total decay width for high- 
energy photons with positive polarization in a lin¬ 
early varying pseudoscalar background is [29] 



(31) 


In order to evaluate an effect of such decays we 
use the same model, as in previous section. We 
assume a layer where the CS vector is pointed 
along the radius. We also assume that the total 
flux of positive polarized photons outgoing from 
the layer is Nq(u>). After photons with positive 
polarization propagate inside the region governed 
by MCS electrodynamics their number should de¬ 
crease as 

N(t, w) = N 0 e (32) 


where 7 = - 1 , stays for the Lorenz factor of 

the particle. In this section we are going to find 
an order of magnitude of the decaying process. 
For the simplicity let us consider photons moving 
along the pseudoscalar gradient. Dispersion law 
for these particles is 


UJ ~\l k+2 + ^ + 2 C ’ 

from eqs. ©• We have neglected here a plasma 
effective photon mass m*. 

Now we will find the 7 factor. The group ve¬ 
locity reads, 


dui Woj 2 — (u 
dk~ u - ±£ < ’ 

(33) 

1 2 w - £ 

7 VT^ C 

(34) 

And we get 


N(t,u) = N 0 e 

(35) 


It is important to remind that there is a threshold 
for the described decay 
„_,2 


u> ~ 


c 


(36) 


This may be a large suppression factor for pho¬ 
tons of positive chirality. Let us assume for the 
sake of discussion that photon energies are of the 
order w ~ GeV and £ is > keV. Then one can eas¬ 
ily see that for t > 10 -9 s most of photons with 
positive chirality have decayed. This time (10” 9 s) 
corresponds to a distance scale ~ 10 cm, which is 
a quite small number for astrophysical object. So 
we have to conclude that if inside the medium 
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there is a strong enough gradient of pseudoscalar 
background, photons with positive polarization 
will decay in e + e - giving increasing the lepton 
pairs going and suppressing almost completely 
the number of photons with positive chirality (of 
course reversing the sign of the gradient the same 
argument applies to the opposite chirality). 

For the process to be possible in a dense fermion 
medium the resulting fermions must have ener¬ 
gies above the respective Fermi sea levels (recall 
that they are slightly different for the two fermion 
‘species’). 


8 Possible implications 


From the previous discussion we have seen that 
there are new potentially interesting effects tak¬ 
ing place in a medium where a pseudoscalar con¬ 
densate with a non-zero gradient is present. Whether 
they are relevant for real physical systems, possi¬ 
bly in compact stars, depends on different circum¬ 
stances. Some of them can be assessed on general 
grounds but others depend on precise astrophys- 
ical details. 

Before entering into the discussion it may be 
interesting to understand the magnitude of the 
relevant parameters b and 

A hypothetic axion condensate related to an 
axion of the Peccei-Quinn type may well couple 
to electrons although its value is largely arbitrary 
but expected to be of order y e x where 

y e is the corresponding Yukawa coupling and f a 
the axion decay constant. a(x ) is the varying ax¬ 
ion field. For neutrons the corresponding coupling 
would be enhanced by a factor ~ 20 4- 30. The 
main uncertainty here is the value to assume for 
(i d x a(x )). This quantity also determines £. 

A pion condensate could make influence on 
electrons with a coupling of 10 -3 less than for 
photons which is related to the tt —> e + e“ decay [30]. 
We follow the paper [21] and refer the reader to it 
for more details. As for neutrons, the influence of 
a varying pion condensate would be substantially 
larger and it can be estimated to be of order 


b 


— 9nNN 


(i d x n(x )) 

U 


(37) 


Thus b is potentially large in such a situation. 

Well known estimates indicate that the Fermi 
energy is of the order of 0.10 MeV for the degener¬ 
ate electrons inside a white dwarf and of the order 
of 50 MeV for the degenerate neutrons in a typical 
neutron star. Clearly in both cases tp «m and 


therefore a non-relativistic approximation such as 
the one used to establish the energy levels is rea¬ 
sonable. 

Possibly the clearest effect predicted is the 
splitting of the Fermi sea in two separate sets 
corresponding to different fermion ‘species’ that 
obey different dispersion relations. In the non- 
relativistic approximation the two Fermi levels 
are separated by the quantity 2b. Under sufficient 
gravitational pressure the two levels will tend to 
equilibrate with a subsequent copious emission of 
photons with characteristic frequency ~ b. De¬ 
pending on specific details these sudden bursts 
could correspond to low frequency radio waves. 
The qualitative estimations of the rate of this in¬ 
teresting process will be done elsewhere. 

The possible consequences due to the change 
in the transmission and reflection properties of 
photons are harder to estimate. In stars it takes 
a long time for radiation to reach surface. For in¬ 
stance, in the Sun the photon diffusion time scale 
is ~ 10 5 years [3Tj. Neutron star are much denser 
(p ~ Po), so the mean free path of photons is 
very small (a full picture about the structure of 
neutron stars can be found in m\) and it is ac¬ 
cepted that, after an initial stage dominated by 
neutrino emission, during most of the life of the 
star heat transfer proceeds through thermal dif¬ 
fusion in the highly degenerated neutron gas and 
radiative cooling is only effective in the outer lay¬ 
ers of the star (crust and atmosphere). The inter¬ 
nal temperature is nearly constant thanks to the 
high thermal conductivity of the neutrons, a con¬ 
sequence of their extremely large mean free paths 
as befits a nearly degenerate Fermi gas. The ac¬ 
tual rate of cooling of a neutron star depends very 
substantially on the opacity of the outer layers 
(‘heat blanketing envelope’ [32] )• Thus it seems 
that the mechanisms discussed in this work could 
be relevant for cooling if large pseudoscalar gradi¬ 
ents are present near the star’s crust, assuming of 
course that the commonly accepted cooling mech¬ 
anisms are correct. 

White dwarfs are constituted by a highly de¬ 
generated relativistic electron gas and the mech¬ 
anism of cooling in their quasi-steady state pro¬ 
ceeds very similarly [33]. Electron-electron inter¬ 
actions are clearly subdominant for same reason 
as neutrons in a neutron star, and transport is 
dominated by electron-ion interactions (ions can 
be treated classically as the temperature is com¬ 
parable to their chemical potential, while T <C 
p e ). Not only that, actually most of the white 
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dwarf heat is stored in the positive ions (recall 
that for the ions classically Cy ~ |lVj, whereas 
for the electron degenerate Fermi gas 
Cy ~ N e T/eF < Cy). Photons in the few keV 
range do not interact easily with electrons due 
to the quantum degeneracy of the latter having 
a Fermi energy of the order of 0.10 — 0.20 MeV 
but they do interact easily with ions and radiative 
cooling is mostly relevant in the external layers of 
the star. 

In any case the mechanisms discussed here 
work in the direction of retarding the process of 
radiative cooling inasmuch as the latter is rele¬ 
vant. Not only a large fraction of the photons 
may be reflected; half of them may actually be 
unstable if the right conditions are given and may 
rapidly decay to fermion-antifermion pairs help¬ 
ing populate the electron sea. 

We have considered other possible influences 
of the distorted dispersion relations, such as changes 
in the fermion conductivity, but they are too small 
to be taken into account. 

9 Conclusions 

In this work we have investigated how the appear¬ 
ance of a pseudoscalar condensate may change the 
properties of photons and fermions inside a stellar 
matter. 

While the formal aspects of the work presented 
here are well founded, the practical relevance of 
the present study hinges on a number of hypoth¬ 
esis. In particular, we have to assume that par¬ 
ity may be spontaneously broken due to the high 
density, to the presence of some axion condensate, 
or both. We have to accept that several domains 
of these characteristics form in the central part of 
the star, if our model pretend to make a contribu¬ 
tion to the compact star inner processes. At the 
very least there should be one domain surrounded 
by an external crust where parity is not broken. 

We also have to assume that the characteristic 
scale of these domains and also the intermediate 
transition regions are much larger than typical 
wave length of photons (and electrons) present 
in the star. This last assumption seems guaran¬ 
teed.The paper considers the situation when the 
pseudoscalar profile is given by two constant val¬ 
ues at x < X- and x > x+ with a linear interpola¬ 
tion in the intermediate region. Certainly it would 
be more realistic to consider a setup where the 
gradient varies smoothly at length scales much 
longer than the photon / fermion wavelength. In 


this case it can be interpolated in piecewise line 
approximation, i.e. the result could be found by 
combining several domains with constant CS vec¬ 
tor in the spirit of Fig. 1 and corresponding con¬ 
volution of photon reflection/transmission coef¬ 
ficients. We postpone this interesting case to a 
forthcoming paper. 

The appearance of a pseudoscalar condensate 
in nuclear matter at high baryon densities has not 
yet been observed but it is predicted theoretically 
in a solid way [9]. On the other hand different au¬ 
thors have speculated with the existence of axion 
condensates (of a Bose-Einstein type) with radii 
ranging from 10 km to 10 4 or 10 5 km. These are 
very relevant sizes for stellar physics. 

Then if these hypothesis hold the mechanism 
of stabilization of the star and even the very pro¬ 
cess of cooling can be affected by the presence of a 
CS vector induced by a varying pseudoscalar con¬ 
densate. There is little or none model dependence 
in the predictions of this phenomenon. The tran¬ 
sition regions are described by Maxwell-Chern- 
Simons electrodynamics and their consequences 
can be worked out independently of the micro¬ 
scopic details of pion/axion condensation. How¬ 
ever in order to get a numerical estimate of the 
modified cooling rate the present mechanism would 
require a knowledge of the distribution of the dif¬ 
ferent domains, at least in average, and more im¬ 
portantly an estimate of the value of the pseu¬ 
doscalar gradients relevant in the present context 
of stellar plasmas. 

Among the effects that we have predicted we 
list the following: (a) Photons and fermions get a 
distorted spectrum. The former have non-trivial 
reflection coefficients by layers where the pseu¬ 
doscalar density, whatever its origin may be, varies. 
Inasmuch as radiative processes are relevant, i.e. 
particularly in regions close to the surface of the 
star if a substantial pseudoscalar gradient is present, 
they can be severely affected retarding the cool¬ 
ing of the star, (b) A similar effect is present for 
fermions helping them to build up larger densi¬ 
ties in the star’s inner regions but its relevance is 
much smaller than for photons. However, one half 
of the fermions are found to be unstable in the 
presence of the axial background and can poten¬ 
tially decay with the emission of a photon. Only a 
small fraction of fermions can actually decay due 
to Pauli’s exclusion principle. This would have 
the opposite effect of helping to cool the star to 
some extent (radiative cooling is inefficient in the 
inner parts of a compact star), (c) A mismatch 
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between the Fermi seas corresponding to the two 
type of fermions may be suddenly rearranged due 
to the star gravitational pressure causing photon 
emission with a well determined spectrum (Fast 
Radio Bursts? j3Jj). Needless to say the practi¬ 
cal relevance of these effects depends mostly on 
the magnitude of the parity breaking parameters 
£ and b. 

Of course if parity is broken inside a star other 
consequences should follow since this would un¬ 
doubtedly modify the equation of state. There 
could also be photon birefringence [28] at the boun¬ 
dary layer and photon instability [29] . All taken 
together could help to detect parity breaking in 
dense (or not so dense) stars. 

If focusing particularly on neutron stars the 
models should describe the cooling rate of these 
objects. In [35] a predicted decay time was ob¬ 
tained from a model of color superconductivity 
and in a recent paper 136] it was shown how the 
crust cooling may depend on the presence and 
properties of nuclear matter. In EH, [28] one may 
find the review of differences between the theoret¬ 
ical predictions and experimental data. All in all, 
it seems fair to conclude that we still do not have 
a good understanding of the cooling of neutron 
stars. 

In the near future a number of new astrophys- 
ical instruments will provide various information 
on neutron stars (see the recent review [39]). In 
particular, the Neutron Star Interior Composition 
Explorer (NICER) [40] to be launched in 2016 is 
expected to discover tens of thousands of neutron 
stars and help understanding their core phenom¬ 
ena and to unravel footprints of parity breaking 
due to pseudoscalar condensation. 
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Appendix A: Solutions of MCS wave 
equation 


We build polarization vectors using the projector 
on the plane transverse vectors k fl , (jj, m 

S\ = 6\D + k"k x ( 2 + C'( x k 2 

-C-k(C X k v + Ck x y, (A.l) 

D = (C ‘ k) 2 — £ 2 k 2 = \S\. 

With a help of the latter equality one can find 
that 

S^e x ^Ck p = De^C a k fi . (A.2) 

Then to our purpose it is convenient to intro¬ 
duce two orthonormal, one-dimensional, Hermi- 
tian projectors 

qpv a 

Ttf = — ± D-J 

= «T = «T; (D > 0). (A.3) 

A couple of chiral polarization vectors for the 
MCS field can be constructed out of constant 
tetrades e„ 

£±{k) = 'n^e tl . (A.4) 

Their properties were thoroughly described in [3T] . 

In order to obtain the normal modes of prop¬ 
agation of the MCS field, let us introduce the ki¬ 
netic 4x4 Hermitian matrix K with elements 

K\u = g\v{k 2 - m 2 ) + is Xva p Ck 0 -, 

K Xv = K* x . (A.5) 


We obtain the general solution for £ • x < 0 from 
the relations (IA.3I) . (IA.4I) 

I<^e±(k) 

5^{k 2 - m 2 t ) + Vd (tt^ - tt^)] e±(fc) 


= ( k 2 — ml ± 


D ) e£(k). 


(A. 6 ) 


Appendix B: Solution of modified Dirac 
equation 


We start from the Dirac equation written in mo¬ 
mentum phase space: 


{YPp - in- = 0 ; 

It is convenient to consider expression: 


7 ° 7 V = 


-« 7 ° 7 1 7 ° 7 1 7 2 7 3 


• 2 3 

= — 1 7 7 = 


oq 0 

0 (7 1 


= -pi ; 


(B.7) 


(B. 8 ) 
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Using this expression and multiplying the first 
equation by 7 0 we get, 

(7°7 /i p /J — 7°m + = 0. (B.9) 

Now we introduce two projectors, 

I ± 5y 


P+ = 


if>± = 'P±ip- 


(B.10) 


One can use these projectors to simplify (IB. 91) 
with a help of the next expressions, 


[0-1,70] = 0; [07,71] =0; 

{07,72} = 0; {0-1,73} = 0, 

namely, 

P+(7°7 M 4V - 7° m + ba^ip = 0 
(po - atipi - 7 °m + 6 ) 0 + 

-a±pj_4>- = 0 ; 

and 

(po - atipi - 7 °m - 6 ) 0 _. 

-a±pj_i/j + = 0; 


(B.ll) 


(B.12) 


(B.13) 


where a* = 7 ° 7 ®, and a±p± = CI 2 P 2 + &3P3- 
To solve our equations (IB.12HB.13l) we multiply 
(IB. 121) by a±p± and get, 

(a±pj_)(p 0 - ct\pi - 7 °m + b ) 


= 


p\ 


>+• (B.14) 


Using this relation one can easily write, 


(po - otipi - 7 °m- b)- 


Pl 


(B.15) 


x (p 0 - aipi - 7 °m + b) - a±p±] 0 + = 0 ; 
and for 0 _, 0 + get equations, 

{pi - p 2 - m 2 - 6 2 

±26(aipi + 7 °m)) 0 ± = 0 (B.16) 

To solve this equation we express 0+ as 0+ = 
( ^ new no ^ a Bon we rewrite (IB. 161) as, 


Oo - P 2 - m 2 

— b 2 ± 2 bm)4>± ± 26p 1 oy£± = 0 
(Po - P - ™ 

— 6 2 =F 26m)£± ± 2bpi<Ji(j)± = 0 
From these equations it is easy to get, 

[(^-P 2 - m 2 -fe 2 ) 2 
—46 2 m 2 — 46 2 p 2 ] 0± = 0, 

wherefrom the dispersion law is, 

Po = P 2 + m2 + 6 2 ± 26^/m 2 + p 2 . 


(B.17) 

(B.18) 

(B.19) 


Thus the Dirac field ip(x) may be written in form, 
0(®) = u A (p)e~ ips+plAXl , (B.21) 


A=± 


where, 
PIT = 


\Jpo — P± + 36 2 — to 2 ± 26^/pp — p\ + 2 b 2 — 2m 2 . 


Appendix C: Transmission/reflection 
through the boundary 

Equation Q in a half-space 27 > 0 enforce v to 
satisfy the following conditions 


(C.22) 


k2k3~\-iuj\/u} 2 — k 2 . „ 

V2 +T± — ui 2 -k% 

^2^3 0J 2 — k 2 , ~ 

— oj 2 — k% V 3~^± 5 

~ ujk^—ik2'\J<jj 2 — k 2 L „ 

— 2(u> 2 -k 2 A ) ’ 

uik 3 Tih 2 ^ ui 2 — k\ „ 

V 0—^± 2(ui 2 — k 2 ^) V 3 —^±, 

V2L-^± , 

. £3 1 ? 3L^ ■ 

Thus we have the solutions in both half-spaces, 
and we should now match them on the boundary. 
If we believe that all contribution to the vector 
field A are continuous, the integration over 27 
from — e to e will give us the next relations m 


~{L) 
UnJ. = 


-rrUo^. + 


Ujk3 


uk2 ~ 
u) 2 — k 2 2- 


oj 2 -k^ U 3- _^ 

<jjks^ik2 \J <jJ 2 — k 2 ± 


~(±) _ k~i_ ~ 

U 0 -+ — ~ 2(uj 2 -k 2 ) M °-» ~ 2 (u 2 -k 2 ± ) 

U)k2±ik3yj uj 2 — k 2 , „ 

2 (uj 2 -k 2 ± ) U2 - > - 


U3- 


~{L) 


U2-> = - 


ko ~ oj ko, ~ 


k-2 k 3 


-X±) 


ui 2 -k 

2 -k% 


VM3- 


COk2 z fik3 yjU) 2 — k' 2 L „ 

~ 2(t o 2 -k 2 ± ) U °- 

k2k3^fioj\/cu 2 — k 2 . „ 

H 2(u, 2 k 2 ± ) - U 3^' 


U 2-> ~ 2{uj 2 — k 2 ^ ) U 2—> + 


~ (L ) /Co ~ aifco ~ 


~(±) u 2 —ki - , uik 3 ±ik 2 Wui 2 -k 2 _ 

U 3-^ — 2{u: 2 -k 2 ± ) U3 ^ _l 2 (oj 2 -fe^) U °- 

. A;2fc3±icd-v/a; 2 — ^ 

^ 2 (w 2 -fe^) M2 -*" 

-,(A) 


where Each component of 

the incoming amplitudes has its own transmission 
coefficient ED, 


^A-t — 


2 fc? 


u 0 1 ;,CS 


0 A) 


(C.23) 


(B.20) 
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Appendix D: Boundary effects for 
electrons 

If we are interested in electrons escaping the parity¬ 
breaking domain, we use the next solutions, 

V>rO) = J2 u A ^{p)e~ m+PlAX1 
A=± 

+ UA^{p)e- ip& ~ PlAXl \ (D.24) 

A=± 

Mx) = w{p)e~ m+PWX1 . (D.25) 


and assume that this particle is on the mass-shell 
Pi = Pia(p)- In this case after the reflection of 
the boundary we would get 


= (0, -foiA -Pio + b2 ) u ^ Q 


2bp 1A ui 


b 2 + (pia + Pio) 2 b 2 + (pia + Pio) 2 


) 


and 


^reflected. = 4b 2 pj A + (p\ A - pj 0 + 6 2 ) 2 

initial {b 2 + {P 1 A + Pio) 2 ) 2 


(D.29) 


Here the first term in ipi stands for the falling 
electrons, the second term is for the reflected ones. 
ip 2 describes electrons penetrating through the 
boundary and propagate in accordance to the Dirac 
equation with pio = \Jp\ — p\ — m 2 . The Dirac 
equation of our system may be used to obtain 
matching conditions, we consider the small area 
near the boundary, 

- m - 7 / %7 5 )V’lxi=^ = 0 (D.26) 

Using our choice of 6-vector, we get, 


(h'd r - 7 1 60 (-x 1 ) 7 5 )^|^ = 0; (D.27) 


Continuity of ip is the first matching condition. 
Since ip is continuous, previous expression leads 
to, 

diip\ Xl=+E - = -*&7 S ^Ui=o- (D.28) 

We now use the (ID. 241 ID.251) forms of solutions 
and rewrite matching conditions in components, 

Pia(wao<- - wao->) - Piow AO = bw A2 ; 

Pia(«ai^ - MA1-S-) - PioWA.1 = bw A3 - 
Pl A (u A2 +- - WA2-!.) - P10WA2 = bw A0 \ 

Pia(u 3 A<- - U A 3_>.) - p w w A3 = bw A1 . 


wao<— — 


UA l<h- = 


U A 2<- = 


2bp 1A u A2 ^ - {p\ A - p \ 0 + 6 2 )mao-» . 
b 2 + (p 1A + Pio) 2 

2bp 1A u A3 ^ - (pi A - Pio + 6 2 )mai-» . 
b 2 + (piA + Pio) 2 

2bp 1A u A0 ^ - (pf A - p \Q + 6 2 )ma 2 -> . 

6 2 + (piA + Pio) 2 
2bpi A u A i^ - (p\ A - p 2 w + b 2 )u A3 


UA ^ b 2 + (pia +Pio) 2 

To give some quantitative effects let us try to 
describe the process of reflection of electrons. We 
deal with wave functions ip, however, only \ip\ 2 
has the physical meaning. Let us for simplicity 
take the propagating electron with 


= (0, ui, 0,0) 
















